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LONG-TERM  GOALS 

The  goal  is  to  investigate,  theoretically  and  through  analyzing  existing  data,  sea  surface  physics  and 
air-sea  exchange  in  extreme  winds.  This  is  a  collaboration  between  Ed  Andreas  and  Kerry  Emanuel. 
Our  underlying  motivations  are  improving  predictions  of  tropical  cyclone  intensity  and  structure  and 
developing  guidelines  for  planning  an  eventual  field  experiment  to  observe  the  air-sea  drag  and 
enthalpy  exchange  in  high  winds.  .  Ultimately  these  goals  require  our  developing  physics-based 
parameterizations  and  theoretical  constraints  for  turbulent  air-sea  fluxes  in  extreme  winds.  One  focus 
will  be  on  the  role  that  sea  spray  plays  in  transferring  heat,  moisture,  and  momentum  across  the  air-sea 
interface  in  high  winds. 

OBJECTIVES 

1 .  Continue  analyzing  data  sets  collected  in  high  winds  (e.g.,  HEXOS,  FASTEX,  CBLAST,  Duck)  to 
deduce  surface  fluxes  and  develop  parameterizations  for  the  air-sea  fluxes  of  enthalpy  and  momentum 
that  begin  to  probe  the  behavior  of  the  sea  surface  in  hurricane-strength  winds. 

2.  Undertake  theoretical  work  to  identify  processes  near  the  air-sea  interface  in  extreme  winds  that 
affect  the  air-sea  exchange  of  enthalpy  and  momentum.  Develop  physical  constraints  for  these 
processes  and  tentative  parameterizations  for  them. 

3.  Do  sensitivity  studies  using  various  ocean  stonn  models  to  evaluate  how  the  parameterizations  for 
air-sea  coupling  that  we  develop  affect  predictions  of  tropical  cyclone  intensity  and  structure  as  well  as 
the  responses  of  the  ocean  and  atmosphere  to  strong  forcing. 

4.  Using  theory  and  models,  combined  with  inferences  about  surface  fluxes  from  the  first  three 
objectives,  quantify  the  sensitivity  of  a  storin’  s  inner  and  outer  structure  and  the  evolution  of  that 
structure  to  assumptions  about  surface  fluxes  and  other  environmental  factors. 

APPROACH 

This  work  is  theoretical  and  analytical;  it  has  no  experimental  component.  Andreas  is  the  only  NWRA 
participant  but  will  be  collaborating  with  Kerry  Emanuel  of  MIT. 
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A  main  emphasis  of  this  work  is  on  how  sea  spray  mediates  the  air-sea  fluxes.  Microphysical  theory 
establishes  how  rapidly  spray  droplets  can  exchange  heat  and  moisture  in  a  given  environment. 

Theory  also  predicts  how  sea  spray  production  should  depend  on  wind  speed  and  how  spray  droplets 
should  be  distributed  in  the  near-surface  air.  The  analytical  part  involves  developing  parameterizations 
for  the  various  spray  transfer  processes  by  simplifying  model  results  or  by  synthesizing  various  data 
sets  and  observations.  Checking  the  parameterizations  against  available  data  is  also  another  aspect  of 
what  I  call  analytical  work. 

As  just  one  example  of  this  recent  theoretical  and  analytical  work,  I  have  developed  a  parameterization 
for  the  salt  flux  to  the  ocean  surface  that  is  mediated  by  sea  spray.  To  my  knowledge,  no  one  has 
recognized  before  that  spray  can  affect  the  ocean’s  buoyancy  by  adding  salt  to  the  surface.  This 
conclusion,  however,  is  a  necessary  extension  of  Andreas  and  Emanuel’s  (2001)  concept  of  re-entrant 
spray.  They  recognized  that  only  spray  droplets  that  are  flung  into  the  air,  cool,  and  then  fall  back  into 
the  sea  can  affect  the  net  enthalpy  flux  across  the  air-sea  interface.  But  similarly,  droplets  that  are 
flung  into  the  air,  evaporate  some  of  their  water,  and  fall  back  into  the  sea  are  saltier  than  the  surface 
water  from  which  they  formed.  Their  evaporation  and  re-entry  into  the  ocean  constitutes  a  surface 
source  of  salt. 

In  my  recently  published  bulk  flux  algorithm  for  high-wind,  spray  conditions  (Andreas  et  al.  2008),  we 
modeled  the  total  air-sea  fluxes  of  latent  (HLT)  and  sensible  (HsT)  heat  as  follows: 


HL  T  +  aQL  , 

(la) 

HsT  =  Hs  +  PQS  -  (a  -  y)QL  . 

(lb) 

Here,  HL  and  Hs  are  the  interfacial  latent  and  sensible  heat  fluxes,  which  we  model  with  the  COARE 
Version  2.6  bulk  interfacial  flux  algorithm  (Fairall  et  al.  1996).  The  QL  and  Qs  are  “nominal”  spray 
fluxes  that  we  compute  with  my  full  microphysical  spray  model  (Andreas  1989)  and  knowledge  of  the 
spray  generation  function  (Andreas  2002)  by  integrating  over  all  droplet  sizes  relevant  to  the  spray 
heat  transfer.  That  is,  QL  and  Qs  include  a  lot  of  microphysical  calculation  that  are  theoretically 
based  but  occur  “off  stage”  in  the  context  of  (1). 

Finally,  the  a,  P,  and  y  are  small,  non-negative  coefficients  that  we  use  to  tune  (1)  to  data.  Andreas  et 
al.  (2008)  used  data  from  HEXOS,  the  study  of  Humidity  Exchange  over  the  Sea,  and  FASTEX,  the 
Fronts  and  Atlantic  Storm-Tracks  Experiment — two  very  good  high-wind  experiments — to  obtain 
a=1.5,  P  =  10.5  ,  and  y=  0.2  .  In  a  modeling  sense,  the  total  fluxes  represented  as  the  left  sides  of  (1) 
would  serve  as  the  lower  flux  boundary  condition  for  an  atmospheric  model. 

The  latent  heat  flux  at  the  air-sea  interface  results  from  evaporation  and,  therefore,  produces  a  salt  flux 
to  the  ocean.  If  the  fractional  surface  salinity  is  s,  the  interfacial  latent  heat  flux  Hl  is  associated  with 
an  interfacial  salt  flux  of 


sHl 

Ml"®)  ’ 


(2) 
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where  Lv  is  the  latent  heat  of  vaporization  of  water. 


Andreas  et  al.  (2008;  also  Andreas  and  DeCosmo  1999,  2002)  compute  QL  in  (la)  from 


Ql  =  PQlOoK  ,  (3) 

where  Qi.(r0)  is  the  spray  latent  heat  flux  contributed  by  all  droplets  of  initial  radius  ro,  and  ri0 
(  =  1.6  pm )  and  rhi  ( =500 pm  )  are  the  smallest  and  largest  droplets  that  contribute  appreciably  to  this 
integral.  Only  the  droplets  that  fall  back  into  the  sea,  however,  contribute  to  the  spray  salt  flux.  These 
are  the  larger  droplets;  suppose  they  have  radii  of  at  least  rcut,  a  cut-off  radius  that  depends  on 
environmental  conditions  such  as  wind  speed  and  ambient  temperature  and  relative  humidity.  Thus, 
using  the  same  a  that  we  obtained  for  (la) — because  I  have  no  other  means  to  evaluate  a — my 
estimate  of  the  spray  salt  flux  is 


F..lw  =  .  “s  'PA PH  •  (4) 

Lv(l  -  s)  Jrcut 

Figure  1  shows  my  calculations  of  the  total  salt  flux,  Fsaitj  from  (2),  where  here  HL  is  the  total 
measured  latent  heat  flux  (Hl.t)  from  the  HEXOS  and  FASTEX  data.  The  figure  also  shows  the  spray 
salt  flux,  FsaitjSp  from  (4).  The  measured  salt  flux  starts  larger  than  the  spray  salt  flux  but,  because  at 
low  wind  speed  it  depends  approximately  linearly  on  wind  speed,  does  not  increase  very  rapidly. 
Although  the  spray  salt  flux  starts  out  low,  because  it  increases  faster  than  the  square  of  the  wind 
speed,  it  will  eventually  equal  and  dominate  the  total  salt  flux  in  stonn  winds.  This  spray  salt  flux  will 
clearly  increase  buoyancy  mixing  in  the  ocean  in  high  winds  and  could  be  an  important  air-sea 
coupling  mechanism  in  storms. 

This  has  been  just  one  example  of  my  approach  in  this  project:  to  use  theory,  data  analysis,  and 
interpretation  to  better  understand  sea  surface  physics  in  very  high  winds. 

WORK  COMPLETED 

The  highlight  of  this  year  was  the  publication  of  my  flux  algorithm  for  high-wind,  spray  conditions, 
Andreas  et  al.  (2008),  which,  as  I  described,  is  based  on  analyses  of  HEXOS  and  FASTEX  data  and 
builds  on  Andreas  and  DeCosmo’s  (1999,  2002)  analyses  of  just  the  HEXOS  data. 

Andreas  et  al.  (2006)  recently  reported  their  analysis  of  the  largest  data  set  ever  used  to  evaluate  the 
von  Karman  constant,  k,  in  the  atmospheric  surface  layer — 553  measurements.  We  concluded  that  k  is 
a  bit  lower  than  the  canonical  value  of  0.40;  our  value  is  0.387 ±0.003  . 

Because  k  appears  prominently  in  all  bulk  flux  algorithms,  I  want  to  update  my  flux  algorithm  to 
reflect  this  new  value.  But  any  such  changes  must  be  made  carefully  because  k  occurs  everywhere  in 
the  mathematics  of  the  atmospheric  boundary  layer;  many  empirical  coefficients  and  functions 
therefore  depend  on  the  value  of  k  assumed  in  the  original  data  analysis.  Before  I  update  my  flux 
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Fig.  1.  Calculations  of  the  total  salt  flux,  Fsalur ,  using  (2)  and  measurements  of  the  total  latent  heat 
flux  ( HLj )  during  HEXOS  and  FAS  TEX,  and  corresponding  estimates  of  the  spray  salt  flux  from 
(4).  The  measured  salt  flux  is  around  10~6  kg  m2  s~2  and  increases  slowly  with  the  neutral-stability, 
10-m  wind  speed,  Umo-  The  spray  salt  flux  starts  at  about  Iff8  kg  m2  s~'  for  Umo  =  5  m  s~2  but 
increase  to  10~6  kg  m  2  s~2  for  Umo  =  20  m  s~2  because  it  increases  faster  than 

the  square  of  the  wind  speed. 
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algorithm,  I  wrote  a  manuscript  that  develops  a  rational  way  to  implement  a  new  value  of  the  von 
Karman  constant  in  many  of  the  equations  used  to  describe  the  atmospheric  boundary  layer  (Andreas 
2008).  That  manuscript  is  in  the  final  stages  of  review. 

Finally,  I  have  been  collaborating  with  Gerrit  de  Leeuw,  Chris  Fairall,  and  several  others  to  prepare  a 
manuscript  that  reviews  our  current  understanding  of  sea  spray  generation  and  dispersion  (de  Leeuw  et 
al.  2008).  The  spray  generation  function,  usually  denoted  dF/dr0 — the  number  of  droplets  produced 
per  square  meter  of  sea  surface  area,  per  second,  per  micrometer  increment  in  the  droplet  radius  at 
formation,  ro — is  crucial  to  all  my  work  (Andreas  2002).  For  example,  it  is  one  of  the  variables  hidden 
in  Ql  and  Qs  in  (1)  and  in  QL(r0)  in  (3)  and  (4). 

Figure  2  shows  a  summary  figure  that  I  have  contributed  to  this  manuscript.  It  depicts  the  spray 
volume  flux,  (4  n r03  / 3)  dF /  dr0 ,  for  a  1 0-m  wind  speed  of  20  m  s”1  for  the  functions  that  Andreas 

(2002)  concluded  were  most  reliable  and  for  several  new  estimates  of  the  spray  generation  function. 
Clearly,  we  still  have  some  work  to  do  to  understand  spray  generation  to  better  than  half  an  order  of 
magnitude. 
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Fig.  2.  Various  estimates  of  the  spray  generation  function  expressed  as  a  volume  flux  fi.e., 

[4  nr*  /  3jdF/dr0]  in  terms  of  the  droplet  radius  at  formation,  r0,  for  a  10-m  wind  speed  (Uw)  of 
20  m  s'.  Andreas’s  (2002)  recommended  functions  were  from  Monahan  et  al.  (1986),  Andreas 
(1992),  Fair  all  et  al.  (1994),  and  Andreas  (1998).  To  these,  I  have  added  functions  from 
Martensson  et  al.  (2003),  for  water  temperatures  of  5°  and  25°C;  and  from  Zhao  et  al.  (2006),  for 
wave  ages  (f$)  of  0.2  and  1.2.  The  data  from  Anguelova  et  al.  (1999)  come  from  a  wind-wave  tunnel 
in  which  the  equivalent  10-m  wind  speed  was  estimated  to  be  20.4  in  s~]  and  predict  higher  fluxes 
than  all  other  expressions.  Finally,  the  SPANDEX  data  (C.  W.  Fair  all,  2008,  personal 
communication)  come  from  another  study  in  a  wind-wave  tunnel  with  an  equivalent  10-m  wind 
speed  of  25  nt  s'.  Most  of  the  functions  agree  that  the  spray  volume  flux  is  between  10~ 12  and 
10  "  m3  m  2  s  1  pm  for  radii  of  10  pm  and  smaller;  most  of  the  functions  also  agree  that  the  volume 
flux  has  a  pronounced  peak  between  10  9  and  10  m3  m  s  1  pm  for  droplets 

with  radii  around  100  pm. 


RESULTS 

To  develop  a  fast  parameterization  for  our  bulk  flux  algorithm,  Andreas  et  al.  (2008)  assumed  that 
spray  droplets  with  initial  radius  near  50  pm  were  good  indicators  of  the  total  spray  latent  heat  flux, 
a  Ql  .  They  therefore  parameterized  this  spray  latent  heat  flux  as 

«Ql  =PwLT 


1  - 


r(^f,5o) 

50  pm 


VL(u.) 


(5) 
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Here,  pw  is  the  density  of  seawater;  VL  is  a  wind  function  that  depends  on  the  friction  velocity,  up  and 
r  is  the  droplet  radius  as  a  function  of  time,  t,  where 

r(t)  =  req  +  (r0  -  req)eXp(-t/Tr)  .  (6) 

In  this,  req  is  the  droplet  radius  at  equilibrium,  a  function  of  temperature,  relative  humidity,  salinity, 
and  initial  droplet  radius,  ro.  Also  in  (6),  t,-  is  the  e-folding  time  for  a  droplet  of  initial  radius  ro  to 
reach  equilibrium.  Finally,  in  (5),  Tf^o  is  the  atmospheric  residence  time  for  droplets  that  started  with 
radius  of  50  pm.  We  base  calculation  of  req  and  xr  on  Andreas’s  (2005)  fast  microphysical  algorithms. 

Andreas  et  al.  (2008)  evaluated  VifuJ  from  their  partitioning  of  the  HEXOS  and  FASTEX  data  into 
interfacial  and  spray  contribution  through  (1).  They  obtained 

VL(u„)  =  1.10xl0-V'22  ,  (7) 

which  gives  Vl  in  m  s_1  for  u*  in  m  s  . 


Because  of  the  similarity  between  (3)  and  (4),  I  can  presume  that  the  spray  salt  flux  follows  scaling  as 
in  (5).  That  is,  I  hypothesize  that,  for  a  fast  flux  calculation,  droplets  of  ro  =  50  pm  also  are 
bellwethers  of  the  spray  salt  flux  and  write 


L  salt,sp 


=  Pv 


vl  -Sy 


1- 


r(^f,5o) 

50  pm 


Vsalt(U0 


(8) 


Thus, 


Ysalt  (U*  )  = 


L  salt,sp 


^  s  ^ 
v  1  s  y 


1- 


r(Tf,5o) 

50  pm 


n3 


(9) 


where  Fsa|Lsp  is  the  spray  salt  flux  shown  in  Fig.  1. 

Figure  3  shows  Vsait(u*)  computed  from  (9)  for  the  HEXOS  and  FASTEX  data.  The  best-fitting  line 
through  the  data  is 


Vsalt(u.)  =  8.01xl0-8u211  ,  (10) 

which  gives  Vsait  in  m  s^1  for  u*  in  m  s-1.  Notice,  (10)  has  almost  the  same  u*  dependence  as  VL  in  (7) 
but  is  slightly  below  that  line.  I  conclude  that  the  droplets  left  out  of  the  Fsa|Usp  computation,  those  for 
which  rlo  <  r0  <rcut ,  do  not  contribute  very  much  to  the  spray  latent  heat  flux,  a  QL . 
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Fig.  3.  The  spray  salt  flux,  FsaitySp,  that  was  obtained  from  the  HEXOS  and  FAS  TEX  data  and 
shown  in  Fig.  1  is  parameterized  here  as  in  (8).  That  is,  this  plot  shows  the  derivation  of  the  wind 
function  Vsalt  (ut)  =  8.0Ix  1 0  8  ill" ,  which  is  the  straight  line  fitted  through  the 

data  on  this  log-log  plot. 


I  will  incorporate  (8)  and  (10)  into  the  next  version  of  my  bulk  flux  algorithm. 

IMPACT/APPLICATIONS 

The  turbulent  air-sea  flux  algorithm  that  I  have  developed  has  three  features  that  are  not  all  present  in 
any  other  air-sea  flux  algorithm:  It  explicitly  recognizes  two  routes  by  which  heat  and  momentum 
cross  the  air-sea  interface,  the  usual  interfacial  route  and  the  spray-mediated  route;  it  has  been  verified 
against  data;  and  it  is  theoretically  based  and,  therefore,  can  be  extrapolated  to  high- wind  conditions. 
My  recognizing  that  evaporating  spray  can  also  add  salt  flux  to  the  ocean  surface  and  my  developing  a 
parameterization  for  this  flux  is  now  a  fourth  feature  that  no  other  air-sea  flux  coupler  has. 

Although  I  have  tested  this  algorithm  against  in  situ  data,  we  still  need  to  see  if  it  improves  predictions 
of  ocean  stonn  structure  and  intensity.  I  am  currently  working  with  Kerry  Emanuel  to  address  this 
issue.  But  I  have  also  had  a  long  collaboration  with  Will  Perrie  and  his  colleagues  at  Bedford  Institute 
of  Oceanography.  In  fact,  they  have  done  the  most  work  in  testing  my  algorithm  in  ocean  storm 
simulations  and  have  documented  their  finding  in  Perrie  et  al.  (2004,  2005,  2006),  Zhang  et  al.  (2006), 
and  Zhang  and  Perrie  (2008).  Generally,  they  find  that  including  the  spray  heat  fluxes  in  their  coupled 
mesoscale  model  gives  better  predictions  for  the  intensity  of  extratropical  storms  when  central  pressure 
and  maximum  surface-level  wind  speed  are  used  as  metrics  of  storm  intensity. 
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TRANSITIONS 


Besides  my  journal  articles  and  conference  presentations  that  describe  my  work  on  air-sea  exchange  in 
high  winds,  I  have  developed  a  software  “kit”  that  contains  the  instructions  and  FORTRAN  programs 
necessary  to  implement  my  bulk  flux  algorithm.  Version  3.1  is  the  current  version  of  that  code. 
Version  3.2,  soon  to  be  released,  will  include  the  new  modules  for  computing  the  interfacial  and  spray 
salt  fluxes.  I  distribute  this  kit  to  anyone  who  asks  about  using  my  algorithm. 

Another  vehicle  for  transitions  is  my  membership  on  the  American  Meteorological  Society’s 
Committee  on  Air-Sea  Interaction.  For  example,  on  behalf  of  our  committee,  during  a  special  session 
on  transitioning  research  to  applications  at  the  2008  AMS  Annual  Meeting  in  January,  I  gave  an 
invited  talk  on  success  stories  and  lessons  learned  in  air-sea  interaction  research.  Secondly,  I  am  the 
co-chairperson  of  the  program  committee  for  the  16th  Conference  on  Air-Sea  Interaction,  to  be  held  in 
January  2009  during  the  Annual  Meeting  of  the  AMS  in  Phoenix.  In  that  role,  I  arranged  for  several 
sessions  relevant  to  the  subject  of  my  current  research  for  ONR.  Namely,  that  conference  will  have 
two  session  (with  12  talks)  on  Sea  Surface  Physics  and  two  sessions  (with  1 1  talks)  on  Tropical  and 
Extratropical  Storms. 

Finally,  I  am  a  member  of  James  Mueller’s  Ph.D.  thesis  advisory  committee.  James  is  Fabrice  Veron’s 
student  in  the  College  of  Marine  and  Earth  Studies  at  the  University  of  Delaware  and  is  working  on  a 
topic  very  close  to  my  own  research — spray  generation  and  dispersion  and  the  role  that  spray  plays  in 
air-sea  heat  and  momentum  exchange.  I  have  been  learning  things  from  James’s  research;  and, 
hopefully,  he  has  found  my  advice  useful. 

RELATED  PROJECTS 

I  am  just  finishing  a  one -year  project  funded  by  the  Mineral  Management  Service.  Kathleen  F.  Jones 
of  the  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory  is  the  PI  and  has  been  funding 
me  as  a  subcontractor.  Our  objective  is  to  develop  guidelines  for  predicting  spray  icing  on  permanent 
platforms  (usually  drilling  platforms)  in  the  waters  around  Alaska.  Spray  icing  is  a  hazard  to  both 
personnel  and  equipment  during  high-wind  events  with  sub-zero  temperatures.  That  is,  the  conditions 
of  interest  in  the  spray  icing  project  overlap  some  of  the  conditions  that  are  important  in  this  ONR 
project.  The  two  projects,  thus,  mutually  leverage  each  other.  For  the  spray  icing  project,  we  have 
been  developing  equations  for  predicting  the  vertical  profile  in  spray  concentration  as  a  function  of 
droplet  radius  from  what  I  know  about  the  sea  spray  generation  function. 
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